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The status of ab initio microscopic calculations of the 2 H(p, 7) 3 He and 3 He(p, e + z/ e ) 4 He 
reactions is reviewed. The methods used to generate accurate nuclear ground- and 
scattering-state wave functions, and to construct realistic electro-weak transition oper- 
ators are described. The uncertainties in the theoretical predictions, particularly those 
relevant to the p 3 He weak capture, are discussed. For the dp radiative capture, the 
theoretical results are compared with the TUNL data in the energy range 0-100 keV. 

1. Introduction 

In the present talk I will review the progress made in the last couple of years in ab initio 
microscopic calculations of the 2 H(p, 7) 3 He [0] and 3 He(p, e + z/ e ) 4 He [0 capture reactions. 
These reactions provide a sensitive testing ground for models of nuclear interactions and 
currents. 

The outline of the talk is as follows. I will first describe the correlated-hyperspherical- 
harmonics method used to obtain accurate wave functions for the A=3 and 4 bound and 
scattering states. Next, I will discuss the model for the nuclear electromagnetic current, 
and will review our current understanding of the radiative capture reaction 2 H(p,7) 3 He. In 
the third part, I will discuss the model for the nuclear weak current, and present results for 
the proton weak capture on 3 He. This process has lately received considerable attention, 
due to the Super-Kamiokande collaboration measurements of the energy spectrum of 
electrons recoiling from scattering with solar neutrinos 0. 

The theoretical description of few-nucleon capture reactions constitutes a challenging 
problem from the standpoint of nuclear few-body theory. For example, since the A=3 
and 4 bound states are approximate eigenstates of the magnetic moment operator, the 
corresponding transition matrix elements between these and the initial dp (or dn) and 
n 3 He states, occurring in the radiative captures, vanish due to orthogonality. As a result, 
these capture processes are extremely sensitive to: (i) small components in the wave 
functions, particularly the D-state admixtures generated by tensor interactions, (ii) many- 
body terms in the electro-weak current operator, and (iii) P-wave capture contributions. 
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2. Wave Functions 

The nuclear Hamiltonian used in the calculations reported here has the form 

# = E^- + E% + E v vk , (i) 

i z " l i<j i<j<k 

where the two-nucleon interaction is the Argonne fig (AV18) model and the three- 
nucleon interaction is the Urbana-IX (UIX) model ||. This Hamiltonian predicts reason- 
ably well the low- lying energy spectra of systems with A < 8 nucleons in "exact" Green's 
function Monte Carlo calculations ||. The experimental binding energies of the trin- 
ucleons and a particle are exactly reproduced, while those of the A=6-8 systems are 
underpredicted by a few percent. This underbinding becomes (relatively) more and more 
severe as the neutron-proton asymmetry increases. An additional failure of the AV18/UIX 
model is the underprediction of spin-orbit splittings in the excitation spectra of these light 
systems. These failures have in fact led to the development of new three-nucleon interac- 
tion models. These developments as well as a discussion of a number of issues regarding 
two-nucleon interactions, such as non- localities, unitary equivalence, etc., can be found in 
Friar's and Pandharipande's || contributions to these proceedings. 

While the nuclear interaction models above are simple to write down, bound- and 
scattering-state solutions for light nuclei have proven to be rather difficult to obtain. 
Intense work in this area and the continuing increase in computational capabilities have 
led, by now, to the development of a number of methods, each with different strengths 
and domains of applicability (for a review and an assessment of them, see Ref. ||). The 
Faddeev-Yakubovsky and quantum Monte Carlo methods are reviewed, respectively, by 
Glockle []n| and Pandharipande [§] in these proceedings. Here, I will briefly discuss the 
correlated-hyperspherical-harmonics (CHH) technique, as implemented by the Pisa group 
for the A=3 |ll] and A=4 |12] nuclei, and summarize a number of results obtained for 



the bound-state properties and low-energy scattering parameters. 

In essence, the CHH method consists in expanding the wave functions on a suitable 
basis, and in determining the expansion coefficients variationally. As an example, the 
trinucleon bound-state wave function is written as 

*=E^ E WW, (2) 

a r cyclic ijk 



where the hyperradius p = yxf + yf (x* and Vj are the Jacobi variables) , and the known 
functions Z a (i; jk) are antisymmetric under the exchange j ^ k and account for the angle 
(xj and fi), spin (sj and Sjk), isospin and tjk), and hyperangle (fa = cos _1 Xj/p) de- 
pendence of channel a. Correlation factors, which account for the strong state-dependent 
correlations induced by the nucleon-nucleon interaction, are included in the functions 
Z a (i;jk). The Rayleigh-Ritz variational principle, (5 u ^f\H — E \^) = 0, is used to de- 
termine the ground-state energy Eq and the functions u a (p). Carrying out the variations 
with respect to the w Q 's leads to a set of coupled second-order differential equations, which 
are then solved by standard numerical techniques. 
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The Nd cluster wave function \fr LSJJz (again, as an example), having incoming orbital 
angular momentum L and channel spin S coupled to total JJ Z , is expressed as 



where the term vanishes in the limit of large intercluster separation, and hence de- 
scribes the system in the region where the particles are close to each other and their 
mutual interactions are large. The term ^\ JJz , instead, describes the system in the 
asymptotic region, and contains the dependence upon the i?-matrix elements, from which 
phase-shifts and mixing angles are obtained. The "core" wave function \&c is expanded 
in the same basis as the bound-state wave function ty, and both the R- matrix elements 
and functions u a (p) occurring in the expansion of \&c are determined by making use of 
the Kohn variational principle. Results for the 3 He and 4 He binding energies, and p 3 He 
scattering lengths, obtained with the CHH method, are listed in Table |l|. 



Table 1 

Binding energies, B 3 and _B 4 , of 3 He and 4 He, and p 3 He singlet and triplet S-wave scatter- 
ing lengths, a s and a t , calculated with the CHH method using the AV18 and AV18/UIX 
Hamiltonian models. 



Model 



B 3 (MeV) 5 4 (MeV) a s (fm) 



a t (fm) 



AV18 

AV18/UIX 



6.93 
7.74 



24.01 
27.89 



12.9 
11.5 



10.0 
9.13 



EXP 



7.72 



28.3 



10.8±2.6 m 



8.1±0.5 
10.2±1.5 ill 







3. The Nuclear Electromagnetic Current 

The nuclear current operator consists of one- and many-body terms that operate on 
the nucleon degrees of freedom: 

j(q) = E j! 1} (q) + E iSf (q) + E j?i(q) , (4) 

i i<j i<j<k 

where q is the momentum transfer, and the one-body operator has the standard 
expression in terms of single-nucleon convection and magnetization currents. The two- 
body current operator has "model-independent" and "model-dependent" components (for 
a review, see Ref. ||). The model- independent terms are obtained from the charge- 
independent part of the AV18, and by construction satisfy current conservation with 
this interaction. The leading operator is the isovector "7r-like" current obtained from the 
isospin-dependent spin-spin and tensor interactions. The latter also generate an isovector 
"p-like "current, while additional model-independent isoscalar and isovector currents arise 
from the central and momentum-dependent interactions. These currents are short-ranged 
and numerically far less important than the 7r-like current. Finally, models for three-body 
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currents have been derived in Ref. ||15|| , however the associated contributions have been 
found to be very small in studies of the magnetic structure of the trinucleons [ |l~5|j . 

The model-dependent currents are purely transverse and therefore cannot be directly 
linked to the underlying two-nucleon interaction. Among them, those associated with 
the A-isobar are the most important ones in the momentum-transfer regime being dis- 
cussed here. These currents are treated within the transition-correlation-operator (TCO) 



scheme |l^,|nj, a scaled-down approach to a full N+A coupled- channel treatment. In 
the TCO scheme, the A degrees of freedom are explicitly included in the nuclear wave 
functions by writing 



N+A 



sl[(i + u, 

i<j 



TR 



(5) 



where \1/ is the purely nucleonic component, S is the symmetrizer and the transition 
correlations U? R are short-range operators, that convert iViV pairs into NA and AA 
pairs. In the results reported here, the \1/ is taken from CHH solutions of the AV18/UIX 
Hamiltonian with nucleons only interactions, while the U?j R is obtained from two-body 
bound and low-energy scattering state solutions of the full N-A coupled-channel problem. 
Both 7iVA and 7AA Mi couplings are considered with their values, fi^NA — 3 n.m. and 
/i 7 AA = 4.35 n.m., obtained from data JTEJ. 



4. The pd Radiative Capture 

There are now available many high-quality data, including differential cross sections, 
vector and tensor analyzing powers, and photon polarization coefficients, on the pd ra- 
diative capture at cm. energies ranging from to 2 MeV |T7|-pO|. These data indicate 
that the reaction proceeds predominantly through S- and P-wave capture. The aim here 
is to verify the extent to which they can be described satisfactorily by a calculation based 
on a realistic Hamiltonian (the AV18/UIX model) and a current operator constructed 
consistently with the two- and three- nucleon interactions ||T|. 

The calculated 5-factor in the cm. energy range 0-2 MeV is compared with data in 
Fig. [TJ, while the predicted angular distributions of the differential cross section a (9), 
vector and tensor analyzing powers A y (9) and T 2 o(9), and photon linear polarization 
coefficient P 1 {9) are compared with the TUNL data below 50 keV from Refs. [[T^HU m 
Fig |2|. The agreement between the full theory, including many-body current contributions, 
and experiment is generally good. However, a closer inspection of the figures reveals the 
presence of significative discrepancies between theory and experiment in the S-factor 
below 40 keV, and in the small angle behavior of a (9) and T2o{9). 

The S-wave capture proceeds mostly through the Mi transitions connecting the doublet 
and quartet pd states to 3 He-the associated reduced matrix elements (RMEs) are denoted 
by m 2 and 7714, respectively. The situation for P-wave capture is more complex, although 
at energies below 50 keV it is dominated by the E\ transitions from the doublet and 
quartet pd states having channel spin 5=1/2, whose RMEs we denote as p 2 and p±. 
The Ei transitions involving the channel spin S = 3/2 states, while smaller, do play an 
important role in T-2o(9). 
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Figure 1. The S"-factor for the 2 H(p,7) 3 He reaction in the cm. energy range 0-2 MeV, 
obtained with the AV18/UIX Hamiltonian model and one-body only (dashed line) or both 
one- and many-body (solid line) currents. In the right panel, only the results of the full 
calculation are shown. 



The TUNL [JL9| and Wisconsin p0| groups have determined the leading Mi and Ex 
RMEs via fits to the measured observables. The results of this fitting procedure are 
compared with the calculated RMEs in Table |2|. The phase of each RME is simply 
related to the elastic pd phase shift ||20|| , which at these low energies is essentially the 
Coulomb phase shift. As can be seen from Table [| the most significant differences between 
theoretical and experimental RMEs are found for \p±\. The theoretical overprediction of 
Pi is the cause of the discrepancies mentioned above in the low-energy (< 50 keV) 5-factor 
and small angle a (9). 

It is interesting to analyze the ratio Tei = \Pi/P2\ 2 - Theory gives Tei — 1, while from 
the fit it results that r^i ~ 0.74 ± 0.04. It is important to stress that the calculation of 
these RMEs is not influenced by uncertainties in the two-body currents, since their values 
are entirely given by the long-wavelength form of the E\ operator (Siegert's theorem), 
which has no spin-dependence (for a thorough discussion of the validity of the long- 
wavelength approximation in E\ transitions, particularly suppressed ones, see Ref. 0). 
It is therefore of interest to examine more closely the origin of the above discrepancy. If 
the interactions between the p and d clusters are switched off, the relation tei — 1 then 
simply follows from angular momentum algebra. Deviations of this ratio from one are 
therefore to be ascribed to differences induced by the interactions in the £7=1/2 doublet 
and quartet wave functions. The AV18 /UIX interactions in these channels do not change 
the ratio above significantly. It should be emphasized that the studies carried out up until 
now ignore, in the continuum states, the effects arising from electromagnetic interactions 
beyond the static Coulomb interaction between protons. It is not clear whether the 
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Figure 2. The energy integrated cross section <j(9)/a (47ra is the total cross section), 
vector analyzing power A y (0), tensor analyzing power T 2 o(8) and photon linear polariza- 
tion coefficient P 7 (9) obtained with the AV18/UIX Hamiltonian model and one-body only 
(dashed line) or both one- and many-body (solid line) currents are compared with the 
experimental results of Ref. [IT? 



inclusion of these long-range interactions, in particular their spin-orbit component, could 
explain the splitting between the pi and p^ RMEs observed at very low energy. Note that 
this discrepancy seems to disappear at 2 MeV 

Table 2 

Magnitudes of the leading Mi and E\ RMEs for pd capture at E p = 40 keV. 



RME 


IA 


FULL 


FIT 


\m 2 \ 


0.172 


0.322 


0.340±0.010 


rri4 


0.174 


0.157 


0.157±0.007 


H 


0.346 


0.371 


0.363±0.014 


M 


0.343 


0.378 


0.312±0.009 



Finally, the doublet rri2 RME is underpredicted by theory at the 5 % level. On the 
other hand, the cross section for nd capture at thermal neutron energy is calculated to 
be 578 fib with the AV18/UIX model, which is 15 % larger than the experimental value 
(508±15) /xb f21 |. Of course, Mi transitions, particularly doublet ones, are significantly 



influenced by many-body current contributions. Indeed, an analysis of the isoscalar (jj,s) 
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and isovector (/iy) magnetic moments of the trinucleons |L5| suggests that the present 
model for the isoscalar two-body currents, constructed from the AV18 spin-orbit and 
quadratic- momentum dependent interactions, tends to overestimate /is by about 5 %. 
The experimental value for fiy, however, is almost perfectly reproduced. The present 
model for two-body isoscalar currents needs to be improved. 



5. The Nuclear Weak Current 

The nuclear weak current and charge operators consist of vector and axial- vector parts, 
with corresponding one- and many-body components. The weak vector current and charge 
are constructed from the corresponding (isovector) electromagnetic terms, in accordance 
with the conserved- vector-current hypothesis, and thus have "model- independent" and 
"model-dependent" components. The former are determined by the interactions, the latter 
include the transverse currents associated with A excitation. 

The leading many-body terms in the axial current, in contrast to the case of the weak 
vector (or electromagnetic) current, are those due to A excitation, which are treated within 
the TCO scheme, discussed above. The axial charge operator includes the long-range 
pion-exchange term required by low-energy theorems and the partially-conserved- 
axial-current relation, as well as the (expected) leading short-range terms constructed 
from the central and spin-orbit components of the nucleon-nucleon interaction, following 
a prescription due to Riska and collaborators 



The largest model dependence is in the weak axial current. The NA axial coupling 
constant g* A is not well known. In the quark-model, it is related to the axial coupling 
constant of the nucleon by the relations g* A = (6\^2/5)gA- This value has often been 
used in the literature in the calculation of A-induced axial current contributions to weak 
transitions. However, given the uncertainties inherent to quark-model predictions, a more 
reliable estimate for g* A is obtained by determining its value phenomenologically. It is well 
established by now |24J] that one-body axial current lead to a ~ 4 % underprediction of the 



measured Gamow- Teller matrix element in tritium /3-decay. This small 4 % discrepancy 
can then be used to determine g* A 0. While this procedure is inherently model dependent, 
its actual model dependence is in fact very weak, as has been shown in Ref. [p4 . 



6. The p 3 He Weak Capture 

The hep capture process is induced by the weak interaction Hamiltonian H\y. After 
partial-wave expansion of the p 3 He scattering state, the transition amplitude is written 
as § 

( 4 He|^b 3 He; Sl , S3 ) = G v £ C^>. r jt(q)|*r?3 J *> , (6) 

LSJJ Z 

where V s and j CT (q) are the lepton and nuclear weak currents, respectively, and C^f^j 
denotes products of Clebsch-Gordan coefficients. The study reported in Ref. @ includes 
all transitions connecting the p 3 He S- and P-wave channels to the 4 He bound state. The 
corresponding wave functions are obtained from realistic Hamiltonians consisting of the 
AV18/UIX and older AV14/UVIII [||] interaction models. 
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The calculated values for the astrophysical S-f actor in the energy range 0-10 keV are 
listed in Table |3|. Inspection of the table shows that: (i) the energy dependence is rather 
weak, the value at 10 keV is only about 4 % larger than that at keV; (ii) the P-wave 
capture states are found to be important, contributing about 40 % of the calculated S- 
factor. However, the contributions from D-wave channels are expected to be very small, 
as explicitly verified in 3 D! capture, (iii) The many-body axial currents associated with 
A excitation play a crucial role in the (dominant) 3 Si capture, where they reduce the 
S-factor by more than a factor of four; thus the destructive interference between the one- 
and many-body current contributions, obtained in Ref. ||16|| , is confirmed in the study of 
Ref. 0, based on more accurate wave functions. The (suppressed) one-body contribution 
comes mostly from transitions involving the D-state components of the 3 He and He 
wave functions, while the many-body contributions are predominantly due to transitions 
connecting the S-state in 3 He to the D-state in 4 He, or viceversa. 

Table 3 

The hep S'-factor, in units of 10 -20 keV b, calculated with CHH wave functions corre- 
sponding to the AV18/UIX Hamiltonian model, at p 3 He cm. energies E—0, 5, and 10 
keV. The rows labelled "one-body" and "full" list the contributions obtained by retaining 
the one-body only and both one- and many-body terms in the nuclear weak current. 
The contributions due the 3 Si channel only and all S- and P-wave channels are listed 
separately. 





E= 


=0 keV 


E= 


=5 keV 


E= 


40 keV 




% 


S+P 


% 


S+P 


% 


S+P 


one-body 
full 


26.4 
6.38 


29.0 
9.64 


25.9 
6.20 


28.7 
9.70 


26.2 
6.36 


29.3 
10.1 



The chief conclusion of Ref. || is that the hep S'-factor is predicted to be ~ 4.5 times 
larger than the value adopted in the standard solar model (SSM) |[26|| . This enhancement, 
while very significant, is smaller than that first suggested in Ref. ||27|| . Even though 
this result is inherently model dependent, it is unlikely that the model dependence is 
large enough to accommodate a drastic increase in the value obtained here. Indeed, 
calculations using Hamiltonians based on the AV18 two-nucleon interaction only and 
the older AV14/UVIII two- and three-nucleon interactions predict zero energy S-factor 
values of 12.1 x 10 -20 keV b and 10.2 x 10~ 20 keV b, respectively. It should be stressed, 
however, that the AV18 model, in contrast to the AV14/UVIII, does not reproduce the 
experimental binding energies and low-energy scattering parameters of the three- and 
four-nucleon systems. The AV14/UVIII prediction is only 6 % larger than the AV18/UIX 
zero-energy result. This 6 % variation should provide a fairly realistic estimate of the 
theoretical uncertainty due to the model dependence. 

The implications of these predictions for the SuperKamiokande (SK) solar neutrino 
spectrum are summarized Fig. |3|. The SK results are presented as the ratio of the measured 
electron spectrum to that expected in the SSM with no neutrino oscillations. Over most of 
the spectrum, this ratio is constant at ~ 0.5. At the highest energies, however, an excess 
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10 11 12 13 14 15 
E [MeV] 



Figure 3. Electron energy spectrum for the ratio between the Super-Kamiokande 825- 
days data and the expectation based on unoscillated 8 B neutrinos [fj6] . 
explanation of the symbols. 



See text for an 



relative to 0.5xSSM is seen (though it has diminished in successive data sets). The SK 
825-day data are shown by the points in Fig. |3] (the error bars denote the combined 
statistical and systematic error). In the figure, the ratio of the hep flux to its value in 
the SSM (based on the hep S-factor prediction of Ref. fl6|]) is denoted by a, defined as 
a = (S ncw / Sssm) x -Pose, where P osc is the /zep-neutrino suppression constant. Presently, 
a = (10.1 x 1(T 20 keV b)/(2.3 x 1(T 20 keV b) = 4.4, if hep neutrino oscillations are 
ignored. The lines in Fig. |3| indicate the effect of various values of a on the ratio of the 
electron spectrum with both 8 B and hep to that with only 8 B (the SSM). In calculating 
this ratio, the 8 B flux in the numerator has been suppressed by 0.47, the best-fit constant 
value for the observed suppression. If the hep neutrinos are suppressed by ~ 0.5, then 
a = 2.2. Two other arbitrary values of a (10 and 20) are shown for comparison. It 
appears that the prediction of Ref. is unable to explain the distortion observed in the 
spectrum at the highest energies. 



7. Outlook 

Improvements in the modeling of two- and three-nucleon interactions and nuclear 
electro-weak currents, and the significant progress made in the last few years in the 
description of bound and continuum wave functions, make it now possible to perform 
first-principle calculations of interesting low-energy reactions involving light nuclei. While 
the extension of the CHH technique to treat systems with A > 6 may prove difficult, vari- 
ational and Green's function Monte Carlo methods should be able to deal with them 
effectively, in particular their pre-capture continuum states. Experimentally known elec- 
tromagnetic and weak transitions of systems in the mass range 4 < A < 9 will provide 
powerful constraints on models of nuclear currents. Work along these lines is being vig- 
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orously pursued. 
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